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ABSTRACT 

We analyze a pair of Suzaku shadowing observations in order to determine the 
X-ray spectrum of the Galaxy's gaseous halo. Our data consist of an observation to- 
ward an absorbing filament in the southern Galactic hemisphere and an observation 
toward an unobscured region adjacent to the filament. We simultaneously fit the spec- 
tra with models having halo, local, and extragalactic components. The intrinsic in- 
tensities of the halo O vii triplet and O Vlll Lyman a emission lines are 9.98lj^ gg LU 
(line unit; photons cm~^ s^"^ sr" -1) and 2.mtoTQ LU, respectively. These results imply 
the existence of hot gas with a temperature of '^10^ '^ K to '^10^ '^ K in the Galactic 
halo. Meanwhile, FUSE O vi observations for the same directions and SPEAR C iv 
observations for a nearby direction indicate the existence of hot halo gas at tempera- 
tures of '^10^''^ K to ~10^'^ K. This collection of data implies that the hot gas in the 
Galactic halo is not isothermal, but its temperature spans a relatively wide range from 
^]^q5.o |.q rv^xO^" K. We therefore construct a differential emission measure (DEM) 
model for the halo's hot gas, consisting of two components. In each, dEM/dlogT is 
assumed to follow a power-law function of the temperature and the gas is assumed to 
be in collisional ionizational equilibrium. The low-temperature component (LTC) of 
the broken power-law DEM model covers the temperature range of lO^-SO _ 106 02 K 
with a slope of 0.30 and the high-temperature component (HTC) covers the temper- 
ature range of lO^''^^ — lO'^ '^^ K with a slope of —2.21. We compare our observations 
with predictions from models for hot gas in the halo. The observed spatial distribu- 
tion of gas with temperatures in the range of our HTC is smoother than that of the 
LTC. We thus suggest that two types of sources contribute to our broken power-law 
model. We find that a simple model in which hot gas accretes onto the Galactic halo 
and cools radiatively cannot explain both the observed UV and X-ray portions of our 
broken power-law model. It can, however, explain the intensity in the Suzaku bandpass 
if the mass infall rate is 1.35 x 10"'^ Mq yr~^ kpc~^. The UV and X-ray intensities and 
our broken power-law model can be well explained by hot gas produced by supernova 
explosions or by supernova remnants supplemented by a smooth source of X-rays. 

Subject headings: Galaxy: general — Galaxy: halo — ISM: general — X-rays: diffuse 
background — X-rays: ISM — ultraviolet: ISMs 
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Introduction 



Not only does hot gas (T > 10^ K) reside in our galaxy's disk, but it resides in the halo. (Here 
we use the X-ray astronomy convention which defines the halo as the region above the majority of 
the Galaxy's H I, thus a bove a heig ht of ^ 150 — 200 pc giv en the parameterization of the H i 



distribution by iFerrierd (Il998al ) and lDickey &: LockmanI (Il990l ). although other conventions would 
call the lower part of this region the thick disk.) Utraviolet and X-ray observations indicate that the 
high-latitude sky is covered by hot gas. Absorption by Galactic O vi ions, tracers of T ~ 3 x 10^ K 
gas, is seen in all of the Far Ultraviolet Spectroscopic Explorer ( FUSE) halo survey spectra of 
sight lines that transit the halo and have high signal-to-noise data (IWakker et al.ll2003l ). 1/4 keV 
X-rays, tracers of T 10^ K gas, are also seen in all directions, but some of these X-rays are 
produced locally, either in the Local Bubble (LB) or in the heliosphere, and by external galaxies. 
After the local and extragalactic co ntribution are subtra cted. X-rays are found to come from most, 
if not all, high latitude directions (jSnowden et al.lll998l ). Such a large covering fraction does not 
require that the hot gas forms a smooth layer. In fact, maps of O vi column density and 1/4 keV 
brightness show a mottled or lumpy distribution. In comparison, maps of 3/4 keV brightness are 
far smoother, with the exception of t he North Polar Spur /Loop I region which is bright in both 
1/4 and 3/4 keV X-rays (see maps in ISnowden et al.lll997l ). 



The height of the hot gas has been found from the O vi column density data. The average den- 
sity of O VI ions falls off exponentially with height above the plane a nd has scale-height s of 4.6 and 
3.2 kpc for northern and southern Galactic hemispheres, respectively (jBowen et al.ll2008l ). Although 
it is not possible to calculate the hot gas scale-height from observations of diffuse X-ray emission, it is 
possible to determine whether or not X-rays are produced beyond clouds of neutral or molecular in- 
terstellar gas. Such analyses, dubbed "shadowing" analysis, find 1/4 keV X-rays originatin g beyond 
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distance from iDanlv et al.l Il993l ) . One analysis of the shadowing filament in the Southern hemi- 
sphere (z = 160 lb 20 pc) reported both the intensity of ultraviolet photons emitted by O vi ions 
and the intensity of 1/4 keV X-rays emitted by hotter gas. Coupling the UV and soft X-ray 
regimes was fruitful; it led to the realization that the radiation by hot gas in the halo accounts 
for a significant fraction of the energy injected into the Galaxy at the Sun's galactocentric radius 
(~ 6 X 10^^ erg s~^ kpc~' ^ vs. an energy inje ction rate of ~ 8 x 10^^ erg s~^ kpc~^ due to supernova 



(SN) and pre-SN winds; IShelton et al.ll2007l l 



Higher energy (0.3 — 2 keV) Suzaku shadowing observations were made for the same south- 
ern hemisphere fil ament as was observed with FUSE and the ROSAT All Sky Survey (RASS). 
Henlev &: SheltonI (2008, hereafter Paper I) processed the raw Suzaku data and extracted both the 
local and halo spectra, but primarily used the data in order to analyze the LB and compare with so- 
lar wind charge exchange (SWCX) contaminated XMM-Newton data for the same directions. Here, 
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we combine the Suzaku spectrum of the halo with ROSAT 1/4 keV and FUSE O vi shadowing 
data in order to create the first long-basehne spectrum of a single region of the Galactic halo. We 
extend our spectrum to 1550 A using C IV data for nearby, but not coincident, pointings taken by 
instruments on the Spectroscopy of Plasma Evolution from Astrophysical Radiation (SPEAR) 
satellite. We compare the X-ray portion of the long-baseline spectrum with model spectra for 
collisional ionizational equilibrium (CIE) and non-equilibrium ionization (NEI) plasmas. We test 
isothermal models as well as models that have two thermal components. Although the X-ray data 
are well fitted by models having two thermal components, the combined UV and X-ray spectrum is 
not. Because the long baseline spectrum requires some EM at lower temperatures, we fit it with a 
differential emission measure (DEM) function. Such a broad spectrum is a powerful tool for testing 
phenomenological models for hot gas in the Galactic halo. We compare our results with three 
models, two of which are phenomenological (accreting gas and supernova remnant (SNR)), while 
the third assumes that the qu antity and temperature of hot gas varies smoothly with height above 
the plane (jYao &: Wana 120071 ). We find that the intensity and spectrum of an SNR are consistent 
with the observations, that the predictions for a simple model of accreted gas under-produce the 
UV intensities relative to the X-ray intensities, and that the geometrical model must be modified 
in order to account for the O vi intensity seen in our direction. 

The observations are described in Section 2 The basic assumptions of the analysis method are 
described in Section 3. Section 4 discusses our tests of the most basic models, those using isothermal 
and two-temperature thermal spectra for plasmas in or approaching CIE. In Section 5, we measure 
the halo's intrinsic O vii and O viii line intensities and use them, together with the intrinsic C iv 
and O VI intensities obtained from SPEAR and FUSE observations, to estimate the halo's EM 
distribution as a function of temperature, T, in the range T ~ 10^ *^ — lO^'*^ K. We follow up this 
preliminary investigation by testing various possible DEM models in Section 6. Our best-fitting 
broken-power-law model is given in Section 6.2. The limitations and the physical implications of 
our modeling are discussed in Section 7 followed by a summary in Section 8. 



2. Observations and Data Reduction 

Shadowing studies require a set of observations, with one observation toward a molecular cloud 
or filament and one nearby. In this study, the on -filament observation was toward a dense knot 



in the filament described bv IPenprase et al.l (119981 ) {I = 278.65°, b = -45.30°). The off-filament 



observation was toward an unobscured direction approximately 2° away {I = 278.71°, b = —47.07°). 
The on- and off-filament Suzaku spectra analyzed here are identical to those described in Paper 
I. In what follows, we just use the data from Suzaku's back-illuminated XISl chip, which is more 
sensitive at lower energies than the three front-illuminated chips. Details of the Suzaku observations 
and data reduction were given in Paper I. Specifically, point sources with 0.2 — 4.5 keV fluxes above 
5 X 10~^^ ergs cm~^ were excluded using the data from the prerelease of the second XMM- 
Newton Serendipitous Source Catalogue. 
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The FU SE O VI intensities (Sheltonll2003l : IShelton et alJl2007l l and the ROSAT R12 and R45 



count rates (jSnowden et alj 119971 ) for the same directions are taken from the existing hterature, 



and the readers are referred to these papers for information on the observations and data reduction. 
The SPEAR C iv intensities for the sight hues near ours were given to us by J. Kregenow (2006, 
private communication) . 



3. Basic Assumptions and Analysis Method 

The total diffuse X-ray emission along most high latitude lines of sight is generally attributed 
to three basic sources, namely the LB, the extragalactic power-law (EPL ) background due to 



unresolved distant activ e galactic nuclei (AGNs), and the Galactic halo (jSnowden et al. 



199S 



Kuntz &: Snowdenll2000l ). Each source is then modeled with one or more components in our multi- 
component model fitting to the observational data. Throughout this paper we adopt a thermal 
plasma component in collisional ionization equilibrium (CIE) for the LB e mission. We use a power 



law with a photon index of 1.46 to model the extragalactic background (IChen et al.l 119971 ). The 
normalization of the power-law model is a free parameter to be determined by the fitting, and we 
obtain ~11 and ~8 photons cm~^ s^^ sr^^ keV^^ at 1 keV for the on- and off-filament lines of 
sight respectively. In this paper, we focus on the hot gas in the Galactic halo, and so we investigate 
various models for the halo component. Since the same models for the LB and EPL components 
are always included in our fitting, for simplicity, we hereafter name the entire multi-component 
model only after the name of the model(s) for the halo component(s). In all of our models the LB 
component is unabsorbed, and the halo and EPL components are subject to absorption. Besides 
the three basic X-ray emitting sources modeled in our spectral analysis, two other sources, i.e., 
SWCX and the X-ray emitting stellar population, are known to produce possible contamination in 
our Suzaku observations. But as we shall argue in Section 7.1, the contamination from these two 
sources is small and can be safely ignored. 

As the Suzaku XISl is not well calibrated below 0.3 keV, we truncate our Suzaku spectra at 
0.3 keV. The final Suzaku spectra we analyze cover the energy range of 0.3-5.5 keV. We also 
omit the data between 1.3-2.3 keV, where the Suzaku spectra are contaminated by instrumental 
emission lines. ROSAT spectra are also available for the same sight lines from the RASS. Although 
the ROSAT spectra are of much lower spectral resolution than the Suzaku spectra, they still make 
a good supplement to the Suzaku spectra because they extend the energy range down to ~ 0.1 keV. 
We therefore include in our analysis the ROSAT R12 data, which cover the ~0.1 — 0.284 keV 
energy band. We fit to the Suzaku and ROSAT R12 (hereafter Suzaku+ROSAT) spectra jointly. 
The higher energy ROSAT bands are not included because these energy bands overlap with the 
Suzaku spectra, which are of much higher signal-to-noise and energy resolution. We use XSPEC to 
generate most of our spectral models and fit them to the observed spectra. Following the argument 
in Paper I, we use the Astrophysical Plasma Emission Code (APEC) vl.3.1 (Smith et al. 2001) 
to simulate thermal models that we fit to the Suzaku spectra and the Raymond & Smith (RS) 
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code (jRavmond &: Smithlll977l : lRavmondlll99ll ) to make thermal models that we fit to the ROSAT 
R12 data. During the fitting, the parameters of the RS model components for the ROSAT data 
are tied to the parameters of the corresponding APEC model components for the Suzaku spectra. 
In all cases we use the phabs model for the absorption (where t ypewriter fon t deno tes XSPEC 
commands). The interstellar medium (ISM) abundance table from lWilms et al.l (120001') is used for 



the c hemical abundances of both the thermal plasma and the absorbing media (see iHenley et al 
20o3). 



Following Paper I, we take A'h = 9.6 x 10^° and 1.9 x 10^° cm ^ (|Henlev et al.l 120071 ) as our 



"standard" values for the on- and off-filament hydrogen column densities throughout the paper. 
Those values were e stimated using the 100 /im intensities, /loO) from the all-sky IRAS maps of 
Schlegel et al.l (119981^ and the /loo-to-A^n conversion relation for the southern Galactic hemisphere 



given m 



Snowden et al.l (|2000l ) . The on- filament column density is consistent with the value derived 
from the color excess of the filament, E(B — V) = 0.17 it 0.05 (jPenprase et al.lll998l ), which yields 
iVn = (9.9 lb 2.9) x lO^'^ cm~^ when scaled using the conversion relation given by iBohlin et al 



(|l978l ). Measurements of A^hi, made from observation s of its 21 cm inten sity, are provided by 
the Leiden- Argentine-Bonn (LAB) Galactic H i Survey (jKalberla et al.ll2005l ). which gives Ahi = 



(5.1ib0.1) X 10 and (2.0ib0.1) x 10 cm~ for the on- and off-filament directions respectively. This 
off- filament Ahi value is similar to the Ah value derived from the /loo intensity, but the difference 
between the on- filament Ah and Ahi measurements is significant. This difference is probably due 
to molecular hydrogen. 

For the uncertaint y in the on-filarnent co lumn density, we use the value derived from the color 
excess measurement of iPenprase et al.l ( 19981 ). yielding Ah = (9.9 it 2.9) x lO^'' cm^^. In order to 
account for the uncertainty in the off-filament column density, we take Ah = 2.0 x 10^" cm~^ as 



the upper limit and Ah = 0.5 x 10 cm as the lower limit, following IShelton et al.l (|2007l ). The 



two models upon which our conclusions are based are tested for various on- and off-filament column 
densities. In the upcoming spectral fits, we find that the uncertainties in the on- and off-filament 
X-ray absorbing column densities have a limited effect on the general results. More details of the 
influence of these uncertainties on the specific models are given in subsections 4.1 and 6.2. 



4. Isothermal and Two- Temperature Halo Models 

4.1. Collisional lonizational Equilibrium Models 

Our investigation of the halo models begins with isothermal and two-temperature models, 
which we assume are either in CIE or in NEI. The results of fitting the model in which the halo is 
assumed to have a single temperature (the isothermal model or IT model) to our Suzaku+ROSAT 
data are shown in the first row of Table [TJ The noted error bars delineate the 90% confidence 
intervals. Hereafter, error bars calculated via XSPEC fits can be taken as 90% confidence intervals. 
The large x^; however, indicates that the isothermal model may be too simple for the hot gas in 
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the Galactic halo along our sight line. Using the RASS data, iKuntz &: SnowdenI (|2000l ) reached the 
same conclusion for the halo's hot gas in general. They found that the emission from the halo's 
hot gas is best described by a two-temperature model (hereafter 2T model). Such a 2T model was 
also adopted in Paper I for the analysis of the Suzaku+ROSAT data, and provided good fits to 
the spectra. A recent analysis of the Chandra observation of a nearby edge-on galaxy (NGC 5775; 



Li et al.l l2008l ) also shows that a 2T model is valid for the halo X-ray emission. We tested the 



2T model for the current project using the Suzaku+ROSAT data, and display the results in the 
fourth row of Table [H The 2T model provides a good fit to our Suzaku+ROSAT data. The 



temperatures we g e t (log Ti = 6.12_qq^ and logT2 = 6.50_qq2) are consistent with those found by 
Kuntz fc SnowdenI (|20od l (logTi = e.OGtoio and logTa = 6A6toll). However, the temp eratures 



we ge t for the hot gas in our own Galactic halo are slightly lower than those found by iLi et al 



([200a) for NGC 5775 (logTi 6.4 and logr2 ~ 6.8). These high temperatures may be connected 
with NGC 5775's higher rate of star formation activity. To test the influence of the uncertainties 
in the on- and off-filament column densities, we vary the on-filament A'h between 7.0 x 10^'^ cm~^ 
an d 12.8 x 10^*^ crn~'^ (co rresponding to the range of values derived from the color excess measured 
by lPenprase et al.l (119981 )). For the off- filament direction, we also test a column density as low as 
0.5 X 10^*^ cm~^ and as high as 2.0 x 10^'^ cm~^. We find that the temperature of the hot halo 
component varies by up to 3%, the temperature of the warm halo component varies by up to 20%, 
and the temperature of the LB component varies by up to 45% from those for the nominal A^h 
values. The derived O vii and O viii intensities (see Section 5.1) vary by up to 22% and 9%, 
respectively. The uncertainties in the on- and off-filament X-ray absorbing column densities have 
little effect on the general conclusions made about the halo's hot gas in this paper. 



4.2. Ionization State of the Halo's Hot Gas 

Here we examine the ionization state of the halo gas by comparing CIE and NEI models for 
the halo emission. As discussed in the previous section, when the ROSAT R12 are included in the 
CIE fitting, we use different CIE models for the different datasets (i.e. RS code for ROSAT and 
APEC for Suzaku). It is, however, impossible to follow this procedure when using NEI model(s) 
because a RS-code-based NEI model is not available in XSPEC. As a result, we fit our models with 
NEI component (s) to the Suzaku spectra only. The parameters of the LB component are fixed at 
the values found from the previous CIE model fitting to the Suzaku+ROSAT data. This is because 
the LB component is mainly constrained by the ROSAT R12 data, which are not included in the 
fitting here. We also fit corresponding CIE models to the Suzaku spectra only (with the LB model 
parameters fixed to the previously determined values; see the second row of Tabled]) and use these 
results for comparison with those of the NEI modeling. 

We first experiment with NEI modeling by testing an isothermal halo. We replace the CIE 
(apec) halo component shown in row 2 of Table [1] with an NEI (nei) component. The NEI model 
has one more parameter r = net, where Ue is the electron density and t is the time since the heating. 
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CIE is reached when r > 10^^ cm~^ s ( Masailll994l ). The fit results are displayed in row 3 of Tabled) 
The best isothermal NEI halo model is similar to the best isothermal CIE halo model in that their 
temperatures are similar, and with a r = s cm the NEI model is approximately 

in CIE. It should be noted that the shape of the NEI halo component is somewhat constrained 
by the fact that the LB component is fixed. However, the NEI model provides a better fit to the 
Suzaku spectra. While examining isothermal halo models may be informative, such models are too 
simplistic. So we progress to two-temperature modeling. 

In preparation for the 2T NEI modeling, we first establish a comparable CIE model, i.e. a 
model that is fit to the Suzaku data alone, but whose LB parameters are taken from the fit to the 
Suzaku and ROSAT data. The parameters for this model are listed in row 5 of Table [1] . We then 
replace one of the CIE halo components with NEI component, while fixing the LB and the other 
halo component parameters to the values listed on row 5. The resulting best fit models are listed 
on row 6 and 7 of Table [H The low value of r in the first halo component in row 6 suggests that 
some of the halo gas may be in the process if ionizing. This is seen again and more strongly when 
we allow both halo components to have NEI (see row 8). 



5. Constraining the EM Distribution of the Halo's Hot Gas using UV and X-ray 

Emission Lines 

5.1. The O VII and O viii X-ray Emission Lines 

The thermal properties of the hot gas in the Galactic halo can be constrained by emission line 
measurements using the halo O vi and C iv intensity measurements obtained from other sources 
and our halo O vii and O viii measurements from the Suzaku data. We will outline an EM 
distribution that spans a temperature range of 2 dex, i.e. T ^ 10^" — lO'''^ K. In Paper I, the O vii 
and O VIII line intensities were measured for both the LB and the Galactic halo. Here, focusing on 
the halo component, we measure the O vii triplet (~570 eV) and O viii Lya doublet (~650 eV) 
line intensities again, using a different method from Paper I. We begin our measurements with 2T 
CIE model that was fit to the Suzaku data, i.e. row 5 of Table [TJ This model provides a good fit to 
our spectra. Although later in this article we show that the 2T CIE model is unable to explain the 
O VI and C iv observations, the accuracy of the measurement of the O vii and O Vlll intensities 
is mainly determined by the goodness of the fitting to the Suzaku data rather than the physical 
meaning of the model. 

The earlier 2T CIE modeling yielded the temperatures and EMs of the hot halo gas. From 
these values we can calculate the intensity of the chosen emission line or complex, /, from 

' = r.^<^^ I « 

where / n^dl is the emission measure, EM, and e(T) is the emission coefficient. The factor of 
nn/ne adjusts for the fact that the values of e(T) tabulated in the APEC database are normalized 
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using n^nH rather than n^. We use the APEC emission coefficient for consistency with our earher 
fitting, in which we used the APEC model. 

The APEC database hsts hne emissivities for a finite number of temperatures. In cases where 
the temperature we are interested in is between two tabulated temperatures, we interpolate to 
obtain the emissivity at our temperature of interest. To calculate the halo's intensity in the O vii 
triplet, we use Equation ([1]) and the temperatures and EMs obtained from the CIE model fitting. 
We include the contributions from the resonance, forbidden, and intercombination lines. In the 
case of the 2T CIE model we sum the contributions of the two halo components. We calculate 
the O VIII emission line intensity in a similar way, including both components of the Lya doublet. 
The O VII and O viii line intensities obtained from fitting the 2T model to the Suzaku spectra are 
shown in the first row of Table [21 

We check our measurements of the O Vli and O vili line intensities using the following in- 
dependent method. We add four 6 functions to our 2T model to represent the O Vii and O viii 
emission: two unabsorbed 5 functions represent the LB oxygen emission and two absorbed 5 func- 
tions represent the halo's oxygen emission. We also "turn off" the oxygen line emission in the 
APEC spectra for the halo components, and fit this new model to our Suzaku data with all of the 
other parameters fixed at their previously determined values. The energies of the two emission 
lines are also free parameters to be determined by the fitting. The best-fitting halo oxygen line 
intensities are shown in the second row of Table El Note that the O Vii and O viii intensities listed 
in Table [2l are intrinsic intensities; in effect, the observed intensities have been deabsorbed with 
respect to the absorp t ion du e to intervening material along the sight line. A similar technique was 



used in iHenlev et al.l (|2007l ) to measure the LB oxygen line intensities. In that case, the oxygen 



emission from the LB APEC model was "turned o f f" by setting the oxygen abundance to zero. 



Here, we refine the technique used in iHenlev et al.l (|2007l ) slightly. Note that the database used 



by the APEC model is made up of two files, apec-vl.3.1_coco.fits and apec-vl.3.1Jine.fits. The 
former is used for calculating the continuum emission and the latter for the emission lines. Before 
running the fitting procedure, we modified the emission line database by setting the emissivities 
for all of the O vii and O viii lines to zero but did not change the continuum database. As a 
result, best-fitting delta functions account for only the oxygen line emission and not the oxygen 
continuum emission. The various techniques for measuring the halo's oxygen line intensities give 
consistent results. Also, as mentioned at the end of Section 4.1, the O Vli and O viii intensities 
are not much affected by the uncertainties in the X-ray absorbing column densities. Henceforth we 
shall use the O vii and O Viii intensities obtained from our standard 2T model. 



5.2. The O vi and C vi UV Emission Lines 



Ultraviolet observations of O vi and C iv resonance line doublet emission from the filament 



region are also available. Using the FUSE observations of the same directions (jSheltonI 120031 ). we 



obtain a de-absorbed intensity from the O vi doublet (AA1032, 1038) of 7750l^QgQ LU (line unit; 
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,20 pjj^ 2_ 



We also 



photons cm~^ s~^ sr~^) with la error bars, assummg an absorbing A^h of 1.9 x 10 
have an off-filament SPEAR observation of the C iv resonance line doublet (AA1548, 1550) to a 
region of size 4' x 4' and centered at a direction (/ = 279.7, b = —47.2) less than 1.0° away from 
our off-filament line of sight (J. Kregenow 2006, private communication). Since we have a C iv 
observation for only one direction, we cannot remove the LB contribution as we do for the other 
lines using shadowing. However, as the LB is measured to have a temperature of ~10^ K, it is not 
expected to emit much in the UV band. This has been confirmed by the O vi observations which 
are sensitive to hot gas of temperature ~3 x 10^ K: the la upper limit of the LB contribution t o 
the O VI doublet is only ~500 LU, or less than 10% of the emission from the halo (|Sheltonll2003l ). 
The ionization potential for C iii^C iv is lower than that of O v— >0 vi, and thus, collisionally 
ionized C iv is sensitive to gas with even lower temperatures. The ~10^ K Local Bubble will 
therefore contribute even less to the C iv emission than it does to the O vi emission, and so we 
attribute all of the C iv doublet emission to the halo. The neutral hydr ogen column density for 
the sight line of the C iv observation is found from the LAB Survey map ([Kalberla et al.ll2005l ) to 
be iVHT = 2.0 X 10^° cm ~^. Using the empirical relat ion NvtT/E(B - V ) = 4.93 x 10^^ cm from 



Diplas &: Savagd ( 19941 ) and the extinction curve from iFitzpatrickl ( 19991 ). we get a deabsorbed C iv 



doublet intensity of 7780 ± 2680 LU from the original observed value of 5790 ± 2000 LU. Since 
C IV emission may also arise from photoionized media, we are going to take this value as an upper 
limit to constrain our model around the temperature of ^ 10^ K. The measurements of the four 
emission features (C iv, O vi, O vii, and O viil) are summarized in Tabled 



5.3. EM Distribution Outlined by the Four Emission Lines 

Here, we calculate the halo's EM as a function of temperature from the intrinsic C iv, O vi, 
O VII and O viii intensities. Equation ([T]), and the theoretical emission coefficient. I n order to 



maintain consistency with our Suzaku+ROSAT modeling, in which the RS d atabase (IRaymond 



1991 



2001 



) was used for photon energies < 0.3 keV and the APEC vl.3.1 database (ISmith et al 
was used for photon energies > 0.3 keV, we take the emission coefficients for the C iv and O vi 
resonance line doublet from the RS database and the coefficients for the O vii triplet and O viii 
Lyo; line from the APEC database. The upper panel of Figure [1] shows the emission coefficients as 
a function of temperature. Because each emission coefficient covers a finite temperature range, we 
cannot simply assign all of the corresponding emission to a single temperature. Instead, for each 
ion, we determine the average emission coefficient for the temperature range for which the emission 
coefficient exceeds 10% of its peak value, e; we multiply this by the range of log T, i.e. A log T, 
then take 47r n^/nH e AlogT to find the average EM per unit logT. The results for the four 
ions, plotted in the lower panel of Figure [H trace out the halo's DEM function. The circles mark 
the temperature at which the emission coefficients peak, the horizontal bars mark the temperature 
range over which the emission coefficient exceeds 10% of its peak value, and the vertical bars are 
error bars calculated from the errors on the intensities. 
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6. DEM Models of the Halo's Hot Gas 



The EM distribution or DEM function outlined by the four ions in Figure [T] show that the 
hot gas in the Galactic halo is not isothermal. The X-ray data also d isallow an iso t herma l halo. 
But, the X-ray data do not disallow a 2T halo. However, as shown in iHenlev et al.l (j2007l ). a 2T 
model fit to the XMM-Newton data for our directions significantly under-predicted (by 3.3o") the 
O VI intensity measured by FUSE. Although the XMM-Newton observations are contaminated by 
SWCX emission, when we repeat the O vi prediction using our Suzaku+ROSAT fit results, we find 
that the discrepancy between the modeled and observed intensities is even larger. Our 2T model 
fit to the Suzaku+ROSAT data predicts an intrinsic O vi intensity of 160 it 27 LU, which is about 
50 times smaller than the in trinsic intensity, 775 0l^nQn LU, which was calculated from the O vi 



intensity observed by FUSE (jShelton et al.l 120071 ) for an assumed Nh of 1.9 x 10^° ^ 



cm 



The fact that the isothermal and 2T models are inconsistent with X-ray and UV measurements 
implies that the isothermal and 2T models may be over simplified. More sophisticated models have 
been proposed by different authors using observations that probe a relative ly wide temp e rature 
range. Based on O Vi, O Vli, and O viii absorption line measurements, lYao &: Wang) (120071 ) 
discussed the non-isothermality of the hot gas in the Galactic halo and proposed a power-law DEM 
model. Assuming exponential temperature and density distributions with respect to the height 
above the Galactic disk, they were able to determine the slope and temperature range of the power 
law. Another power-law DEM mod el covering a temperature range of ~10^'*^ — 10^'^ K has been 
constructed bv IShelton et al.l (|2007l ) for the Galactic hot gas based on FUSE O vi and ROSAT 
R12 (1/4 keV) and R67 (1.5 keV) observations for our on-filament and off-filament directions. As 
these DEM models were proved successful for hot halo gas within certain temperature ranges, we 
test various DEM models to see if they are consistent with our set of observations covering a wide 
range of temperature as shown in Figure [TJ 



6.1. Power-Law DEM Models 



Shelton et al.l (|2007l ) took the differential path length to be a power-law function of temperature 



dl = BT^dlnT for Ti < T < T2. For an isobaric gas, this results in a power-law DEM model of 
the form 

dEM{T) ^ J yj^^ 

otherwise, 



dlogT 



oc < 



(2) 



where a = f3 — 2. IShelton et al.l (|2007l ) fixed the low-temperature cutoff at Ti = 10^ *^ K because gas 
of lower temperature makes negligible contributions to the O vi intensity. For Nh = 2.0 
they obtained a = -0.05 ± 0.17 and T2 = 10^-^ K. 



10^0 cm-2 



The model is also shown in the lower panel of Figure [2] as the dotted line. We can see a 
good agreement between the model and our O vi emission line data point. The agreement is 
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expected, because both the model and the O vi data point are produced using the same O vi 
measurements (the smah discrepancy is mainly due to the d ifference in the oxyg en abundances, 
emission coefficients and absorption column density assumed in lShelton et al.l (|2007l ) and this paper; 
see below for more details). 



In this section we would like to rework the IShelton et alJ model, to see if it is also consistent 
with the Suzaku observations which cover a higher energy band. But before doing this, we would 
like to take the opportunity to improve upon some of their approximations. For example, they 
used very modern a. bundances in calculating the O vi intensity (i.e. 0/H = 4.57 x 10"'* from 
Asplund et al.l (|200J)), while the RS co de they used for the model ing of the ROSAT spectra relied 
upon an older set of abundances from lAnders &: Grevessd (|l989l ). Here we adopt a single set of 
modern and consistent abundances for the modeliii g of al l of the data. Following the argument 
in Paper I, the abundance table from IWilms et al.l (120001 ) is used in this paper with an oxygen 
abundance of O/H = 4.90 x 10~^. In order to predict the O vi and R12 intensities for comparison 
with the measurements, we need the O vi emission rate coefficient, royi{T), and the R12 emission 
rate coefficient, ri2{T), which are equivalent to l/47r nn/ne e(T) in Equation ([T|). We extract these 
from the RS code because both the O vi emission line and the R12 band are at energies lower than 
0.3 keV. The O vi emission rate coefficient, tq vi(^)) is consistent with the e(T) used in the previous 
section for the evaluation of the O vi doublet intensity. Our R12 emission rate coefficient, ri2{T), 
is obtained by convolving the ROSAT R12 response function with spectra of various temperature 
plasmas that were calculated using the RS code. Our R12 emission coefficient, ri2(T), pertains to 
the intrinsic intensit y and we use the deabsorbed R12 count rate in our calculation. In contrast, 
Shelton et al.l (j2007l ) used an emission coefficient, ri2{Nn,T), that accounted for the absorption of 
material along the sight line and used the absorbed halo R12 count rate in their calculation of the 
DEM. Our ri2(T) is equivalent to their R12 emission coefficient, ri2{N}i,T), when A^h is set to 
zero. 

To measure the deabsorbed R12 intensity from the Galactic halo's hot gas, we fit our 2T CIE 
model simultaneously to the on- and off-filament ROSAT spectra extracted from the RASS database 
(jSnowden et al.lll997l ). Here, the RS model is used for all of the LB and halo components because 
we focus on the R12 band. Throughout this paper, the off-filament hydrogen column density is 
taken to be Ah = 1.9 x lO^'' cm~^ except for the testing of the influence of t he uncertainty i n thi s 



value, so readers are reminded to compare our results with those obtained by lShelton et al.l (120071 ) 



for the most similar Ah value (2.0 x 10 cm~ ). The intrinsic R12 emission from the Galactic 
halo's hot gas is then easily calculated using the flt results. Our best-fitting spectrum yields 

Taking the deabsorbed O yi intensity as /q vr ~ 



R12= (3740 ± 450) x 10"^ counts s'^ arcmin'^ 



7750 



+950 



1090 LU, yields Iq vi/R12 = 2.07; 



-0.36 , 
-0.38 • 



10^- 



LU 



Following IShelton et al.1 (|2007|) 



we determine the index of the power-law model using their Equation (8) which we reproduce here: 



^12 j;^:^^n2{T)T-dlnT ■ 



(3) 
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However, we replace their R12 emission rate coefficient, ro vi (-^H ) 7") , with the deabsorbed emission 
coefficient roYi{T) so as to be co nsistent with our de absorbed measurement of the R12 count rate. 



and we use a = /? — 2. Following Shelton et alj (j2007l ). we set the low-temperature cutoff, Ti to be 



IQ^.o then test this model for various high-temperature cutoffs to see if the model is also 

consistent with our Suzaku+ROSAT data. For each high-temperature cutoff, T2, we calculate the 
slope a using Equation ([3]). Since both the value of /qvi and R12 (not just the ratio) are known to 
us, we can also determine the constant of proportionality in Equation ([2|). With this constant, the 
model predictions for the O vii and O viii line intensities are then calculated using the analogs of 
the numerator of Equation ^ . The vq vii and tq viii emission rate coefficients are extracted from 
the APEC database as these two lines have photon energies > 0.3 keV. 



Here, our technique diverges from that of IShelton et al.l . We try several plausible values for 



the high-temperature cutoff (logr2=6.06, 6.24, 6.36, and 6.54). For each value, we determine a 
from Equation ([3]), determine the constant of proportionality for Equation ([2]), and then calculate 
the O VII and O viii line intensities predicted by the power-law model. Table H] lists the a and the 
O VII and O viii intensities for each examined T2 values. The observationally determined intrinsic 
O VII intensity (9.98l}'gg LU) is best modeled by the second case, that having T2 = lO^-^'' K, 
a = 0.54, and /qvii = 10.4 LU. This case slightly overpredicts the O Vll intensity but is within the 
observational error bars. The other cases over- or underpredict the intensity by > 50%. However, 
the model significantly under-predicts the O viii intensity (0.49 vs. 2.66^q3q LU). If we increase T2 
to 10^-^^ K to improve the agreement between the model predicted O viii intensity (2.97 LU) and 
the observed intensity, the model then more severely over-predicts the O vii intensity (18.4 LU). 
The poor correspondence between model and observation can be seen more directly in the upper 
panel of Figure El where the power-law model with T2 = 10^'^^ K is shown as the dotted line, 
together with the four emission line data points. Since the model is derived from the FUSE O vi 
observations, it does match the O vi emission line data point well. The O vii and O viii line 
data points, however, obviously drop away from the model. They and the O VI data cannot 
simultaneously be explained by a single power-law model with any choice of the high-temperature 
cutoff, as demonstrated in Tabled! 



The power-law DEM model of lShelton et al.l (120071 ) was constructed mainly based on the FUSE 



O VI and ROSAT R12 measurements, but we have found that it is impossible to extend the model 
to a higher temperature and make it consistent with the O vii and O viii measurements. We now 
test the power-law model in another way. We ffist fit a power-law model to the Suzaku+ROSAT 
data and then determine if it is consistent with the C iv and O vi measurements. The exponent, 
a, and the high-temperature cutoff, T2, are free parameters of the model to be determined by the 
fitting. The low-temperature cutoff, Ti, is not well determined by the fitting and is therefore fixed. 
We test models with various choices of Ti and present the results in Table [5l While this power-law 
model fits the Suzaku+ROSAT data quite well for all of the choices of the low-temperature cutoff, 
Ti, none of the values of Ti results in predicted C iv and O vi intensities that are both consistent 
with the measurements. For a low-temperature cutoff of Ti = 10^'^^ K the predicted intrinsic C iv 
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intensity is consistent with the observational value, but the predicted intrinsic O vi intensity is 
far too large (46,500 vs. 7750+?^^o LU), whereas for Ti = lO^-^^ K, even though the predicted 
O VI intensity is consistent with the observational value, the predicted C iv intensity is too small 
(181 LU) unless essentially all of the observed C iv emission (7780±2680 LU) is due to photoionized 
gas in high-pressure photoionized regions around hot stars. When the low-temperature cutoff is 
Ti = 10^'^^ K, both the O vi and C iv intensities predicted by the model fall below the observed 
values. Again, this is shown directly in the upper panel of Figure El The dot-dashed line, which 
represents the power-law model that best fits the Suzaku+ROSAT data and has a low-temperature 
cutoff Ti = 10^'™ K, is consistent with the O vil and O viii emission line data points, but is not 
consistent with the O vi and C iv data points in the lower temperature range. 



6.2. A Broken Power-Law DEM Model 

Motivated by the partial successes of our power-law model that was fit to the Suzaku+ROSAT 
data and the power -law model that was found for the O vi and 1/4 keV emission (patterned after 



Shelton et al 



20071 ). we investigate a broken power-law DEM model of the form 




dlogT f Ti\ if To < T < Ts 



Like the power-law model of lShelton et al.l (120071 ). the slope of the low-temperature (Ti < T < 



T2) portion of the broken power-law DEM, ai, is constrained by the FUSE O vi intensity and some 
of the ROSAT R12 count rate, using Equation We fix Ti to the value of 10^'® K (rather than 
Iq5-0 j^-j^ gQ tj^at the low-temperature portion of the model fully covers the temperature regime 
probed by the C iv emission line observation. The slope of the high-temperature {T2 < T < T3) 
portion, 02, is constrained by fitting to the Suzaku+ROSAT data. While for a reasonable break 
temperature (T2 > 10^'^ K) nearly all of the O vi emission is produced by the low-temperature por- 
tion of the DEM, the R12 emission is produced by both portions. Hence, when using Equation ([3|) 
to calculate ai, the denominator on the left-hand side should not be the total R12 count rate, but 
instead the fraction of the total R12 count rate produced by the low-temperature portion. We 
calculate this fraction (and ultimately the shape of the broken power-law DEM) using the following 
iterative procedure. Because we do not know the break temperature, T2, a priori, we repeat the 
procedure using several different values of T2 between 10^'^^ K and 10^'^'' K. We first fit the broken 
power-law DEM model to the Suzaku+ROSAT data with ai fixed at some initial estimate, and 02 
and the normalizations free to vary. From these fit results we calculate the R12 count rates from 
the two portions of the DEM, and insert the R12 count rate due to the low-temperature portion 
into Equation ([3]) in order to calculate a new estimate of ai. We then re-fit our model to the 
Suzaku+ROSATspectra with ai fixed at the new value, and with 02 and the normalizations free 
to vary. We repeat this procedure until the new R12 count rate due to the low-temperature portion 
of the broken power-law model differs from the old one by less than 10%. This method turns out 
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to be quite efficient and stable. The slopes converge within 4 or 5 iterations for a wide range of 
initial values of ai . 

For the Suzaku+ROSAT data, the fit results with the broken power-law DEM model are 
summarized in Table [H The break temperature of the model is well constrained, as it is easy 
to understand that too high or too low a break temperature will make the broken power-law 
model essentially fail for the same reason that the single power-law models failed. Various break 
temperatures around 10^''' K are tested. Besides the x^/dof value, the model predicted C iv 
intensity is used as a second constraint. Of the models we tried, the one with a break temperature 
of lO^ '^'^ K and an a2 of -2.21 is preferred. This model has the smallest /dof as well as the best 
agreement with the observed C iv intensity. 

We also test the effect of the uncertainties in the on- and off-filament X-ray absorbing column 
densities on our BPL model. Using the same uncertainty ranges as those used in Section 4.1, we 
find that the O vii and O Viii intensities differ by 9% and 8% at most, respectively, and the slopes 
of the high-temperature component (HTC) we get are consistently equal to —2.2. The uncertainty 
in the off-filament column density does affect the slope of the low-temperature component (LTC) 
of the BPL and the 1/4 keV intensity predicted by the model could vary by up to ~20%. Both of 
these effects, however, are not significant enough to affect the main conclusions of this paper. 

7. Discussion 

7.1. Possible X-ray Contamination in Our Suzaku Spectra 

The assumption that our Suzaku observations are not severely contaminated by SWCX X-rays 
was supported in Paper I, where it was shown that the foreground oxygen intensities measured 
from the Suzaku spectra are consistent with zero. In addition, the on-filament and off-filament 
observations were completed within a couple of days of each other during a minimum in the solar 
activity cycle. 

Note that what we have called the LB component in our models really accounts for all of the 
foreground emission, including the emission from the LB and that from SWCX, if there is any. 
Thus the O vii and O viii intensities of our LB component (found by fitting our composite model 
to the data) actually provide upper limits on the corresponding intensities due to SWCX. As shown 
in Tables [Tl[5l and[6l all of the models (except for the IT model, which we have already shown does 
not provide a good fit to the spectra) yield similar results for the so-called LB component, with 
T ~ lo^'^^ K and EM ~ 7 x 10~^ cm~® pc. These fit results for the LB component then predict 
foreground O vii and O viii intensities of ~ 0.16 and ~ 0.06 LU, consistent with the results found 
in Paper I using a different method and 1.0 ± 1.1 LU for O VII and O VIII respectively). 

Also, the foreground O vii and O viii intensities predicted by our LB component are much smaller 
than the corresponding intensities due to the hot gas in Galactic halo (~ lO.Ol?;^ and - 2.71JJ LU 
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for O VII and O viii respectively). As a result, even if the small foreground O vii and O viii 
intensities found by our modeling are all due to the SWCX, they should not affect our analysis of 
the hot gas in Galactic halo. 

We now discuss the possible X-ray contamination from the stellar population. In processing 
the raw Suzaku data, we removed sources with 0.2 — 4.5 keV fluxes above a critical flux fc = 
5 X 10~^^ ergs cm~^ s~^. Unresolved X-ray emission from fainter stellar sources with fluxes < 
5 X 10~^^ ergs cm~^ s~^, however, could be mixed with the diffuse X-ray emission from hot gas. 
We therefore estimate the fraction of the observed X-ray emission that is due to stellar sources with 
fluxes less than the critical flux fc- To do this, we used the X-ray luminosity function, for the entire 
stellar population in the solar neighborhood, given in lSazonoy et al.l (12006). To be conservative, 
we use the space density profile for stars with My > 3.5 from lojha et al.l ( 1996 ). Assuming that 
the luminosity function is independent of height, the total number of stars within our Suzaku field 
of view is ^ 

N= [ n{z,L)dVdL= [ n(z)^(L) dLdz, (5) 

J J sm'^ |6| 

where n{z) is the density profile as a function of height z, and ^{L) = dN/dL^ is the X-ray 
luminosity function, Q = 17.8' x 17.8' is the field of view of our Suzaku observations, and b ~ —45° 
is the Galactic latitude of our observations. Since the flux associated with each star is / 

Lsin^ 161 



L 



the total X-ray flux from all of the stellar sources within the field of view is 
^n{z,L)dVdL = 

Jz=0 



ft 



^=°° 1 n 

— n(z)$(L)L— 

L=0 



:dLdz. 



And"^ ' ' " 7^=0 Jl=o 47r'"'~'^' ' sin 161" 
The total X-ray flux from fainter stellar sources with fluxes less than the critical flux is 

1 

=0 Jl=o 

where = 47r/c[2;/ sin 



fs 



L=Lc 



4tt 



n{z)'^{L)Lsm\b\ndLdz. 



(6) 



(7) 



The space density profile for stars with > 3.5 from lOjha et alJ (119961 ) is 

z z 

n{z) oc e ''1 + 0.074e ''2 , 



(8) 



: 260 pc and scale height of the thick disk is /i2 = 760 pc. 
2.4 keV band derived from the RASS for stars in the 



where the scale height of the thin disk is hi 
The X-ray luminosity function in the 0.1 
solar neighborhood is shown in Figure 5 in lSazonov et al.l (|2006l ). Following the form used fro the 
2 — 10 keV X-ray luminosity function given in the same paper, we parameterize the X-ray luminosity 
function for the 0.1 — 2.4 keV band as 



dN 



dlog Lx 



K 



{Lb/L 



X} 



• Lx < Lb 



[{Lb/Lxr\Lx> Lb. 
From Figure 5 in ISazonov et al.l (|2006l ). we estimate Lb ~ 2.0 x 10^^ 



(9) 



-,-1 



02 1.02. We find that 



ergs s ai ~ 0.6 and 
of the stellar X-ray emission was excluded by our point source 
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removal, and only ~ 20% of the stellar X-ray emission remained to be mixed with the diffuse X-ray 
emission from the hot gas. Note that fainter stars have a more extended distribution than the 
brighter ones. Because we use the density profile for stars with > 3.5, we overestimate the 
fraction of stars located at a larger distance from the plane and thus overestimate the fraction of 
stellar X-ray emission that was not removed by our flux cutoff. Also, not taking into account X-ray 
absorption makes the estimated fraction of stellar X-ray emission that is unresolved higher than 
the true value, as absorption tends to increase with distance. Since the total X-ray emission from 
the stellar populatio n was estimated to be comparable with that from the hot gas in our Galaxy 
(jSazonov et al.ll2006l ). and our conservative estimate shows that we have removed at least ~ 80% 
of the stellar X-ray emission, we argue that our X-ray observations of the Galactic hot gas are not 
badly contaminated by the X-ray emission from the stellar population. 



7.2. Comparing BPL DEM Models with 2T Models 

Using all of the FUSE, ROSAT and Suzaku observations available for our sight lines, we 
have successfully constructed a broken power-law DEM model covering a temperature range of 
lO^'^" — 10^ '^^ K. Before we discuss the physical implications of this model for the hot gas in the 
Galactic halo, we point out that the 2T model, the single power-law model, and the broken power- 
law DEM model all give similar values of when fit to the Suzaku+ROSAT spectra, although 
the 2T and single power law are inconsistent with the UV observations. This means that the 
Suzaku+ROSAT data alone are insufficient to distinguish between these various halo models. We 
can further demonstrate this fact with fake spectra generated from one of our better-fitting broken 
power-law DEM and 2T models. We first generate fake on- and off- filament Suzaku and ROSAT 
spectra from the broken power-law DEM model that has a break temperature of lO^ ''^ K (row 3 
of Tabled]) with the XSPEC f akeit command. When we fit the resulting spectra with our broken 
power-law DEM model, we obtain x^/dof = 537.1/535, whereas when we fit them with our 2T 
model we obtain /dof = 538.9/535. We also generate fake spectra from our 2T model. These 
spectra give /dof = 564.7/535 when fit with the T2 = IQ^-^^^ K broken power-law DEM model and 
X^ /dof = 558.3/535 when fit with the 2T model. The fact that both models give similar quality fits 
to both sets of fake spectra shows that the Suzaku spectra cannot distinguish between the broken 
power-law DEM and 2T models. We thus conclude that even though the 2T model provides a good 
fit to the Suzaku+ROSAT data, it may not necessarily be the real physical condition of the hot 
gas in the Galactic halo. As for the single power-law model, which is plotted in the upper panel 
of Figure [2] as a dashed line, we can see that it mimics the high-temperature part of the broken 
power-law model for the temperature range ~ 10^-^ - lO''-^ K covered by the Suzaku+ROSAT data, 
but overproduces the O vi intensity. 
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7.3. Chemical Abundance Used in the Modehng 



Following iHenlev et al.l (j2007l ). we have adopted the interstellar abundance table from lWilms et al 



pood ) for both e mitting and a 



dance table from 



Wilms et al 



those from Anders &: Grevesse ( 



jsorb ing gas in our spectral modeling. However, besides the abun- 



( 2000 ), five othe r abundance tables are ava i lable ii i XSPEC, namely 



19891 ) ■ lFeldmanl(jl992l ^. lAnders fc Ebiharal (jl982l ). lGrevesse fc Sauval 



(| 19981 ) ■ and|Loddery (l2003|). Here we test the T2 = 10^ °^ K BPL's sensitivity to the choice of abun 



dance tables. For these tests, we first chose abundance tables which will be used for both the LB 
and halo components. We then set the low-temperature limit of the first part to be Ti = lO'^'^ K 
and the break temperature to be T2 = lO^ '^^ K. We use the technique described in §6.3 to determine 
the slope of the LTC. All of the other parameters are determined by the fitting, which we repeat 
until the fit results stabilize. The results are shown in Table [71 The gross structure of the BPL 
remains almost the same, regardless of the abundance table used. There are some differences in the 
slope of the second part of the model, but given the error bars, the slopes are relatively consistent 
with each other. We thus conclude that our model results are fairly independent of the choice of 
abundance table, and our discussion of the properties of the hot gas in the Galactic halo based on 
the models is not affected by the uncertainty in the abundances of the emitting and absorbing gas. 



7.4. Two-Component Scenario for the Halo's Hot Gas 

As clearly demonstrated by the spatial differences between the RASS maps in the R12, R45, 
and R67 bands, the 1/4 keV surface brightnes s is very patchy while the 3/4 keV and 1.5 keV maps 



are much more uniform (jSnowden et al.l 119971 ). This difference between the 1/4 keV and 3/4 keV 
emission cannot be ascribed entirely to the heavier absorption in the 1/4 keV band. Instead, 
the halo's 1/4 keV and 3/4 keV X-rays may be produced by different components. The hot gas 
which produces the majority of the halo's 1/4 keV surface brightness is not uniformly distributed 
while the hotter gas that produces most of the 3/4 keV emission is more smoothly distributed. A 
comprehensive discussion of a two-component model, of course, requires observations of multiple 
sight lines, but here we would like to show that our study of one part of the sky is consistent with 
a two-component model. 

Our BPL model can be naturally divided into two parts, with each part being a power law. 
We name these two parts the LTC and the HTC in accord with the temperature ranges they cover. 
The R12, R45, and R67 surface brightnesses made by the different components of the BPL model 
are listed in Table El Before we focus our attention on the R12 and R45 bands, we note that the 
total R67 count rate is dominated by the contribution from the EPL component. This component 
is expected to be fairly isotropic across the sky, in good agreement with the smoothness of the R67 
band RASS map with the exception of the known extra radiation originating from the direction of 
the center of the Galaxy. The emission from the HTC accounts for most of the total R45 count rate 
(i.e., ~ 55%), while the EPL component accounts for the remainder (~ 40%). The latter alone. 
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however, cannot entirely explain the smoothness of the RASS R45 map; therefore, we expect the 
halo's 3/4 keV X-ray emitting gas (which is modeled in this paper as the HTC) to be fairly smooth. 
In the R12 band, the emission from the LTC and the LB make a significant fraction (~ 90%) of the 
total R1 2 count rate. The LT C account for ~ 50% of the non-local emission, which according to 
maps in lSnowden et al.l (jl997l ). has a patchy distribution. Since the R12 count rate from the HTC 



should be fairly c onstant, then th e LTC must be responsible for the patchy appearance of the non- 
local R12 map in lSnowden et al.l (|l997l ). In a simplified picture of the halo, our BPL accounts for 
the two-component nature of the halo, with the hot gas with a patchy distribution being modeled 
by the LTC of our BPL and the hotter gas of much more uniform distribution being modeled by 
the HTC of our BPL. The possible origins of the two kinds of hot halo gas will be discussed in 
Section 7.6. 



7.5. Comparison with the Other DEM Models for the Halo's Hot Gas 



DEM functions covering a similar temperature range have also been constructed by variou s 
authors for the Galactic hot gas (e.g., Ishelton et aD boO?! : Ivao fc Wane! boO?! : lYao et al.l[2009l ). 
As already mentioned in Section 6.1 and Section 6 .2, ou r BPL model is actually motivated by 
the fact that the power-law model of IShelton et al.l (120071 ) cannot be simply extended to a higher 
temperature range and still be consistent with our Suzaku measurements of the O vii and O viii 
intensities. For consistency with our Suzaku observations, a HTC with a slope of must 
be added to a LTC. The result is a BPL model that is consistent with all of the UV and X-ray 
observations. Our BPL model then mainly differs from the power-law model of lShelton et al.l (|2007l ) 
by the addition of a HTC which covers a temperature range of 106-°2-l07-02 K, and accounts for the 
new Suzaku observations which are sensitive to the emission of the hot gas within that temperature 
range. 



Yao &: Wang (120071 ) also derived a power-law DEM function for the Galactic halo's hot gas, 
based on the assumption that the z distribution of both the temperature, T, and density, n, follow 
exponential functions characterized by scale heights of and Setting 7 = hx/hn, their DEM 
function can be written as 

dEM(r) ^ 



dlog(T) 



oc 



(10) 



Yao et al.l (|2009l ) also used this functional form. Their power-law DEM models, constrained by 
the O VI and soft X-ray absorption line observati ons and diffuse soft X-ray emission observations 
for the hot gas i n the directions toward Mrk 421 (jYao fc Wana 120071 ) (dashed hue) and LMC X-3 
(|Yao et al.ll2009l ) (dot-dashed line), are also shown in the lower panel of Figure [2l 



Two significant differences between the power-law models of lYao &: Wangj (j2007l ) and I Yao et al 



([20091) and our BPL model can be noted. First, their power-law models have slopes of 1.2 and 1.0, 
respectively. Not only are they inconsistent with our slope of ~— 2.2 for the HTC of our BPL, but 
they are also much higher than the slope of ~0.3 for the LTC. Also, their power-law models have a 
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much lower dEM(T)/dlog T value in the low-temperature range than our BPL model. However, we 
have used an O vi intensity measurement and ROSAT R12 data to constrain our model, while they 
did not include an O vi intensity in their modeling. Their use of the O vi absorption observations as 
a constraint ensures that th eir model is consi stent with the O vi column density, but not necessary 
the O VI intensity. In fact, lYao et al.l (|2009l ) did note that the O vi intensity predicted by their 
model is much less than typically observed for the halo. They suggested that most of the observed 
O VI intensity could be due to a second phenomenon, one that is not modeled by their calculation: 
a transition temperature region between hot and cool gas. 



The values of dEM{T)/d\ogT near T = lO^-^ K in the models of lYao fc Wand (|2007|) and 



Yao et al.l (1200911 do agree wit h our BPL model in the high-temperature range, as both the power- 



law model of I Yao et al.l (120091 ) and our BPL were constrained using the O vii and O Vlll intensity 
measurements in the respective directions. Also, their use of n{z) and T{z) suggests that the 
modeled hot gas is smoothly distributed (this is further supported by their similar results for the 
two different directions). This is consistent with our interpretation of the T > 10^ '^ K gas, as we 
also suggest that a portion of this gas component may be smoothly distributed. It would be useful 
if additional data sets could be obtained, so that the analysis of the thermal the spatial properties 
of the Galactic halo's hot gas could be repeated for a larger number of directions. 



7.6. Implications of Our Model to the Origin and Distribution of the Halo's Hot 

Gas 

It is widely believed that the hot gas in the halos of galaxies such as our own is due to 
either accretion of the intergalactic medium (IGM) or stellar feedback. The thermal and spatial 
properties of the hot gas resulting from these two mechanisms could be quite different. Therefore, 
some indications of the origin of the hot gas may be obtained by comparing the thermal properties 
deduced from the observations with those predicted by the theoretical models or simulations. In 
addition, given our two-component scenario for the Galactic halo's hot gas, the morphology of each 
component may also give some useful clues regarding the origin of the component. In this section 
we consider the possible origin(s) of the Galactic halo's hot gas by comparing our DEM model and 
emission line intensities with those predicted by a simple accretion model and by SNR simulations. 
The morphology of the X-ray emission expected from these two models is also compared with the 
RASS maps as another indicator of the origin of the Galactic halo's hot gas. 



7.6.1. Accretion model 

Firstly, we consider a simple cooling model for IGM gas that is heated as it accretes onto 
the Galaxy. For simplicity, we assume that the gas thermalizes as it falls though the Galaxy's 
gravitational potential, comes to rest at some distance from the Galactic center, and then begins 
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to cool radiatively. Furthermore, if we ignore subsequent raising, falling, expansion, or contraction 
of the gas parcels, then each parcel's potential energy can be taken to be constant because no 
additional work will be done on the parcel. We estimate the temperature of the gas before it begins 
to cool radiatively, Tq, to be somewhat less than 10^'^ K, assuming that it was heated due to falling 
though the Galaxy's gravitational potential to a Galacto-centric radius of ~8.5 kpc and assuming 
that the electrons equilibrate with the stripped hydrogen and helium ions. Gas of this temperature 
contains O vii and O viii ions and emits 3/4 keV X-rays, whose spectrum is appropriate for 
comparison with observations. 

The internal energy per unit volume of a parcel of gas at temperature T is Uy = ^nkT, where 
n is the number density of particles. The gas is assumed to be fully ionized, so n = ng + n^, where 
rig and Ui are the number densities of electrons and ions, respectively. The gas is assumed to be 
fully ionized. The parcel of gas loses energy at a rate of dU^/dt = — nenjA(T), where A(T) is the 
cooling function. From the internal energy and the loss rate equation, we can determine that the 
plasma's temperature changes at a rate of: 

dT ^ -2nen,A(T) 

dt 3nk ' ^ ^ 



Equation (jlip can be used to determine the EM function of the hot gas that accreted onto the 
Galaxy's halo if we assume that the gas accretes at a steady rate and then begins to cool down 
from the same initial temperature, Tq. We take the accretion rate, dN/dt, to be constant, where 
N is the number of accreted particles within a cross sectional area, A. Suppose that the thickness 
of the accreted layer is / and that the hot gas is produced and cools down in a way that both the 
density, n and the cross sectional area remain constant. This is the isochoric case. In this case, 
dN = nAdl during the time interval dt. If the accretion proceeds for a time span that is longer than 
the cooling time for the temperature regime of interest (T ~ lO'^'® — 10®'^ K), then the gas reaches 
a steady state with respect to temperature. Accreted gas enters the system at a temperature of 
Tq, cools over time, and eventually leaves the temperature regime of interest. However, as a given 
segment of gas is cooling, newer gas replaces it. Thus at any given time, the quantity of gas of any 
given temperature (within the temperature regime of interest) remains constant. The quantity of 
material within a temperature interval, dT, can therefore be calculated from: 



dN 

It 



dN dt 
'dTdT 



dN 3nk 

(12) 

dt 2neniA{T)' ^ ' 



The EM associated with a given interval is dEM = nj,dl, where dl is simply dN/nA. Thus, 
dEM(r) _^2dl _ fne\ /3A;\ / 1 dN\ 1 ^^^^ 



dT ""dT \ni) \ 2 ) \A dt J A{T)' 



This equation can easily be compared with our plotted DEM function given that d log T = 
log{e)dT/T, A(T) for isochoric, solar abundance, CIE gas is tabulated in [Sutherland &: Dopita 



- 21 - 



(jl993l ). and noting that Ue/ni^ 3A;/2, and {1/ A)dN / dt are constants, although, admittedly, {1/ A)dN / dt 



is of unknown value. 

A scaled version of the DEM function for this simple steady state cooling scenario is plotted 
in Figure El together with our broken power-law model. The two curves bear little resemblance to 
each other. Allowing the accr eted gas to have subsqlar a bundances by using the A(T) curve for 



0.1 solar metalicity CIE gas in [Sutherland &: Dopital (|l993l ) would slightly change the curve's slope 



between T = 10^-° - 10^ ° K, but as shown in Fi gure O would not bring the theoretical DEM into 
agreement with the observationally determined DEM. 

While the DEM functions predicted by the simple isochoric accretion model and the approx- 
imation of CIE do not match the BPL DEM function we got from fitting both the UV and soft 
X-ray data of our sight lines, we note that, within the temperature range 10^'" — 10^'^ K, the 



accret ion r nodel pred i cted D EM functions mimic the power-law DEM functions of lYao &: Wane 



(j2007l ) and lYao et al.l (|2009l ) , which successfully explain the soft X-ray emission seen on their lines 
of sight. This suggests that even though the accretion model fails to be a single explanation for 
all of the UV and soft X-ray emissive hot gas in the Galactic halo, it still might be a phenomeno- 
logical explanation for the soft X-ray emissive portion of the hot gas. To test this idea, we fit 
the accretion model DEM to our Suzaku data only. The technique mentioned in Section 6.2 is 
used again to construct a tabulated accretion DEM model. This time, the spectra for 50 differ- 
ent temperatures obtained from the APEC database are weighted by the accretion model DEM 
function (Equation ()13p ). We calculate the accretion DEM spectra this way for a grid of values 
of high-temperature cutoff with the low-temperature cutoff being set to be 10^'° K. The high- 
temperature cutoff and the scaling of the model are then two free parameters to be determined 
by the fitting. Fitting the accretion model to the Suzaku data results in a reasonably good fit, 
with x^/dof = 574.8/533. The accretion rate can also be estimated from the best-fit normalization 
value. Our result is M/A = {U / 23) {mn / A) {dN/dt) = 1.35 x lO'^ Mq yr"^ kpc"^ The factor 
14/23 comes from assuming that the accreting gas is fully ionized with H:He=10:l. If intergalactic 
material falls evenly onto our galaxy across the whole disk, then A = IttB? , where i? = 15 kpc is 
the radius of the disk, and the factor of 2 is for the two sides of the disk. For the whole galaxy 
we then have M = 1.9 Mq yr^^. Technically, this is an upper limit, because some of the X-rays 
observed by Suzaku may come from other sources. 



7.6.2. SNR model 



Next, we consider the possibility that the gas was heated by an explosive event, such as an 
SN explosion. For this c o mpar ison, we use the results of SNR simulations from the series of papers 
by ISheltonI (jl998l . Il999l . 120061 ). The simulations employed a Lagrangian mesh hydrocode with 
algorithms tha t model shock dynamics, NEI and recombination, nonthermal pressure and thermal 
conduction. In IsheltonI ( 2006 ). modeled SNRs located between 130 to 1800 pc a bove the Galacti c 
midplane. The density of the ambient medium at these heights was taken from [j^rrierel ( 1998b ). 
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As we can see in Figure IH both the O Vl and X-ray emission of an SNR are significant only when 
the SNR is younger than ~ 10^ yr. The number of SNRs of age ^ 10^ yr that reside on a typical 
sight line can be estimated from , typical radius (< 100 pc) of the SNR and the SN explosion rate 



(Equation (10) in ISheltonI (|2006l )). Because the typical number is small (< 0.05), the probability 
of encountering two or more SNRs of age ^ 10^ yr on a sight line is tiny and we thus only compare 
our observationally results with those predicted for a single SNR. 

The predicted O vi/R12 ratio for an SNR residing at z ~ 1300 pc (ambient density of no=0.01 
atoms cm~^) has alre ady been compared with the FUSE and ROSAT observations for our directions 



(jShelton et al.l l2007l . although different units and conversions were used in that paper), and the 
conclusion drawn that the observed ratio best matched that of a remnant at an age before the SNR 
formed a dense shell. Here we extend the work by comparing with a greater number of simulated 
SNRs and extending the comparison to the 3/4 keV X-ray band. We examine remnants evolving 
in ambient densities of no=0.2, 0.1, 0.05, 0.02, 0.01, and 0.005 atoms cm~^, corresponding to 
heights of 2:=190, 310, 480, 850, 1300, and 1800 pc, respectively, SN explosion energies of Eq = 0.5 
and 1.0 X 10^^ ergs, and ambient nonthermal pressures of Pnth = 1800 and 7200 K cm~^. The 
integrated O vi, R12, and R45 intensities predicted by various model SNRs are shown in Figured] 
as curves of different colors and line types. Although our FUSE and Suzaku observations may 
cover only a small portion of an evolved SNR, the exact positions of our sight lines relative to the 
possible SNR(s) are unknown. However, we find that the variation of the intensities from one sight 
line to another is fairly small, and the integrated SNR intensities are thus used in comparing the 
SNR models with our observations. Our observationally derived intrinsic O vi intensity, 1/4 keV 
count rate and 3/4 keV count rate are 7750l?^^o LU (from Section 5.2), 3740 ± 450 x 10"^ R12 
counts s~^ arcmin"^ (from Section 6.1), and 92^3^ x 10~^ R45 counts s~^ arcmin"^ (from the 
prediction of our BPL model). Interestingly, the observationally derived O vi intensity and R12 and 
R45 count rates simultaneously match the predictions for the SNR that has no=0.01 atoms cm~^, 
£"0 = 1.0 X 10^^ ergs, and Pnth = 7200 K cm~^ (corresponding to the dashed blue curve in Figured]) 
at an age of ~ 10^ yr. Although the other models are able to match the halo O vi intensity and 
R12 count rate, they are not able to match all observations simultaneously. 

The consistency between the observations and the predictions for a single SNR, however, is 
not what we expected from the discussion in Section 6.3, where we argued that most of the homo- 
geneously distributed 3/4 keV-emitting gas may have different source than the inhomogeneously 
distributed 1/4 keV-emitting gas. Like the 1/4 keV-emitting gas, gas traced by O vi ions is also 
inhomogeneously distributed. Because the number of SNRs encountered on a sight line is small, 
the sporadic SNe should produce a patchy distribution of hot gas. Therefore, it is logical to relate 
the LTC of our BPL model, which produces almost all of the O vi emission and the majority of the 
1/4 keV emission, to the SNR model, and to relate the HTC of our BPL model, which produces 
most of the 3/4 keV emission, to the more uniformly distributed hot gas. For this reason, we 
now compare the O vi, R12, and R45 intensities derived from the LTC with the predictions of the 
simulated SNRs. As shown in Figure H] the intensities derived from the LTC are best matched 
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by the predictions of the SNRs with no=0.02 atoms cm at an age of ~ 1.8 X 10^ yr ( green 
curves). Although the simulated SNRs overpredicts the LTC's R45 intensity, the predicted value 
is still smaller than the halo's total intrinsic R45 intensity and thus leaves room for R45 emission 
produced by more smoothly distributed hot gas. 

While both the O vi and R12 portions of our halo emission observations can be explained by a 
single SNR, the collection of O vi column density and R12 count rate m easurements for the high 
latitude sky do not show a constant ratio between O vi and soft X-rays ([Savage et al.ll2003l ). We 
note, however, that the variation in attenuation from one line of sight to another is not taken into 
account in these surveys nor is the contribution from the local region. Also, SNRs create both O vi 
and soft X-rays, but the ratio of O vi to R12 changes as the remnant evolves. For the simulations 
we discussed, the ratio of O vi intensity to R12 count rate varies by a factor of >100 with age. For 
these reasons, SNRs should not be expected to result in a constant O vi to R12 ratio on all line 
sights across the high latitude sky. 

We would like to end the discussion of SNRs as a possible component of the halo's hot gas with 
two interesting points. First, based on the filament's IRA S 12 nm to 100 um int ensity ratio, IRAS 
60 fxm to 100 fim intensity ratio, and Ca ii kinematics, iPenprase et al.l (|l998l l determined that 
the filament has been heated, probably by a shock, s uggesting that t he fil ament could be a very 
old SNR. Second, we note that the halo R12 maps of lSnowden et al.l (jl998l ) clearly show that the 
filament and our off-filament observation overlap with a roughly circular region of angular radius 
~ 5° that has an elevated R12 count rate and an unusually low R2/R1 ratio. These characteristics 
imply the existence of an unusual feature in this direction. For comparison, the angular radius of the 
simulated SNR can also be estimated from the distance and the predicted radius of the simulated 
SNR. We find that for the 1.8 x 10^ yr old SNR that has no = 0.02 cm'^, Eq = 1.0 x 10^^ ergs and 
Pnth = 7200 K cm~^, the angular radius is ~ 4.0°, which is consistent with the angular size of the 
distribution in the R12 maps. 



8. Summary and Conclusions 

In this paper we analyze the Suzaku spectra of the ISM obtained from observations pointing 
toward and to the side of an absorbing filament at high southern Galactic latitude. We take a joint 
analysis of these data and FUSE and ROSAT observations of the same sight lines, in order to 
constrain the thermal and spatial properties of the hot gas in the Galactic halo. Our main findings 
are as follows: 

1. O VII and O viii emission line features due to the Galactic halo's hot gas are firmly detected 
using our Suzaku shadowing observations. Their intrinsic intensities are 9.98l{:^[| and 2.66l[J:i^5 LU, 
respectively. These observations, together with the FUSE observations of emission from O vi in 
the Galactic halo along our off-filament line of sight (intrinsic O vi intensity = 7750^*^^090 ^"^^ 
SPEAR observations of emission from C iv along a direction less than 1.0° away from our off- 
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filament line of sight (intrinsic C IV intensity = 7780 it 2680 LU), sample hot gas with temperatures 
ranging from ~ lO^ *^ to ~ 10^'^ K. These observations indicate a non- isothermal distribution of the 
hot gas in the Galactic halo, which is coi isistent with the finding o f earlier auth ors who mode l ed the 



halo with two thermal co mponents (e.g. iKuntz &: Snowdenll200d ). and that of lYao &: Wand (j2007l ) 



and lShelton et al.l (|2007l ) who modeled the halo with power-law DEM functions. 



2. We construct a DEM model for the halo's hot gas, in which dEM{T)/dlogT follows a 
broken power-law function of T. This model is consistent with the SPEAR, FUSE, ROSAT, and 
Suzaku observations. The LTC of the broken power law covers the temperature range from lO^ '' K 
to lO^''^^ K, has an index of ai = 0.30 and is mainly constrained by the FUSE O vi intensity 
and ROSAT R12 count rate. The low-temperature cutoff of this component can be extended to 
10^'^ K and is consistent with the SPEAR C IV intensity. The HTC covers the temperature range 
from lO^ ''^ K to lO^'*^^ K, has an index of 02 = —2.21 and is mainly constrained by the Suzaku 
X-ray spectra. If we take the X-ray emission in our sight line to be representative of that of the 
entire halo, then we can estimate the 0.2 — 2.0 keV soft X-ray luminosity of our Galaxy's halo 
to be 3.0 X 10^^ erg s^^. Considering the spatial differences between the 1/4 and 3/4 keV RASS 
maps, we propose that the Galactic halo's hot gas is composed of two components. The higher 
temperature, more uniformly distributed component is represented by the HTC of our BPL, while 
the lower temperature, less uniformly distributed component is represented by the LTC of our BPL. 
Confirmation of this hypothesis, of course, requires observations in more sight lines. 



3. We coinpare o ur BPL DEM model with the power-law DEM models of lYao Wang ()2007l ) 



and lYao et al.l (|2009l ) for the Galactic halo's hot gas toward the directions of Mrk 421 and LMC 
X-3, respectively. Their models assumed an exponential disk scenario. Comparing their power-law 
models with our broken power-law model, we find that the curves are inconsistent. We propose 
that the fundamental reason for this inconsistency is that we included O vi emission information 
in our analysis, while they did not. 

4. We compare our results with the following scenarios: 

(a) A simple IGM accretion and cooling scenario. In this scenario, we assume that intergalactic 
gas is accreted onto the Galactic halo at a constant rate. The hot gas then cools radiatively 
such that a line of sight though the accretion layer samples a range of temperatures. We derive 
the DEM function for the accretion layer, finding that its shape is inconsistent with the broken 
power-law model derived from the observations (see Figure [3]). It under-predicts the UV intensity 
relative to the X-ray intensity. However, with a high-temperature cutoff value of 10^'^ K, the X- 
ray emission predicted by the accretion model is loosely consistent with our Suzaku observations. 
Attributing all of the Suzaku-hmiA X-ray emission to the accretion model yields an accretion rate 
of 1.35 X 10~^ M0 yr~^ kpc~^ or 1.9 Mq yr~^ for the whole galaxy. This is an upper limit because 
some of the X-rays seen by Suzaku may have come from other sources. 

(b) A SNR scenario. We use existing simulations of SNRs evolving at various heights above 
the disk. We find that the observed O VI, 1/4 keV, and 3/4 keV intensities match the predictions 
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of a ~10^ yr old SNR located at a height of 1300 pc above the disk. The predicted angular size 
of such a remnant is consistent with a bright spot on the ROSAT 1/4 keV map in the direction 
of our observations. Because we suspect that a more smoothly distributed source supplements the 
3/4 keV intensity of sporadic explosive events, such as SNRs, we also consider dimmer SNRs. A 
slightly older SNR located nearer to the galactic plane can make the O vi and 1/4 keV photons and 
some of the 3/4 keV photons while 'leaving room for' a smoother source of 3/4 keV photons (see 
Figure E]). The IGM accreted onto our galaxy could be a possible origin of the smoothly distributed 
component of the hot gas. 

We thank J. Kregenow for providing us with the SPEAR C iv intensity measurements for sight 
lines near ours. We would like thank the referee for very helpful comments. This project was sup- 
ported by NASA through the Long Term Space Astrophysics Program under grant NNG04GD78G, 
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Newton Guest Observer Program under grant NNG04GB08G. This research has made use of data 
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Fig. 1.— Upper panel: The C iv (AA1548, 1550), O vi (AA1032, 1038), O vii (triplet at ~0.57 keV), 
and O VIII (~0.65 keV) emission coefficients as a function of the gas temperature. The values for 
C IV and O vi have been scaled down by a factor of 1000 for clarity. The C iv and O vi line 
coefficients are from the RS database, those for O Vll and O Vlll are from the APEC vl.3.1 
database. Lower panel: Galactic halo EM distribution as outlined by the four emission lines. The 
four circles mark the temperature at which the emissivity peaks and the EM per unit log T needed 
to produce the halo's C IV, O VI, O VII, and O VIIl intensities. The horizontal error bars show 
the measurement of the intensities of the lines, and cover the temperature ranges for which the 
theoretical emission coefficients are more than 1/10 of their peak values. The vertical error bars 
are derived from errors on the intensities. 
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Fig. 2. — Upper panel: Two power-law DEM mo dels tested in t his p aper for the hot gas in the 
Galactic halo: a power-law model similar to that of lShelton et al.l (|2007l ) (dotted line) and a power- 
law model fitted to the Suzaku+ROSAT data only (dashed line). Lower panel: O ur best-fitting 



broke n power-law DEM model (so lid lines) in compa rison with the power-law of IShelton et al 



(|2007l ) and the power-law model of lYao Wangl (j2007l ) towards the direction of Mrk 421 (dashed 
line) and that of lYao et al.l ([20091) to wards the direc tion of LMC X -3 (dot-dashed line). (Note that 
we have reduced th e DEM models of lYao &: Wang) and lYao et al.l using the oxygen abundance of 



Wilms et al.l (|200d ).) The C iv, O vi, O vii, and O viii data points from the lower panel of Figure 



[T]are shown in both panels for comparison. 
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Fig. 3. — Comparison between our BPL and theoretical cooling models: Because the {l/A)dN/dt 
coefficient is not known, we set it arbitrarily to a value that allows the accretion DEMs to overlap 
our observationally determined BPL DEM. The scaled DEM functions predicted by simple cooling 
models (dotted and dashed curves) do not resemble our BPL model above a temperature of lO^^^ K, 

showing that our model is not consistent with a simple cooling picture, in which the hot gas is first 
heated to an X-ray emitting temperature of ~ 10®'^ K (see the text for more details) and then 
cools radiatively. The dotted curve is the DEM function predicted by isochoric cooling of solar 
metallicity CIE gas and the dashed one is for isochoric cooling of 1/10 solar metallicity CIE gas. 
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SNR Age(yr) 

Fig. 4. — Observationally derived halo O vi, R12, and R45 intensities (horizontal solid lines in each 
panel, with error bars noted by set of solid lines at the top right of the panel) are compared with 
the predictions for various simulated SNRs (colored curves). The O vi intensities are in unit of 
photons cm~^ s~^ sr~^ and the R12 and R45 intensities are in units of 10~^ counts s~^ arcmin"^. 
The ambient densities (no=0.2, 0.1, 0.05, 0.02, 0.01, and 0.005 atoms cm^^) of the simulated 
SNRs are keyed to the curve colors using a rainbow sequence (red, orange, yellow, green, blue, 
and purple). The line types of the curves (solid, dotted, and dashed) distinguish the SN explosion 
energy (£^0=0. 5, 0.5, and 1.0 x 10^^ ergs) and ambient nonthermal pressure (P„(/t=1800, 7200, 
and 7200 K cm~^). The total halo O vi, R12 and R45 intensities are fairly consistent with the 
predictions of the SNR with no=0.01 atoms cm~^, Eq = 1.0 x 10^^ ergs, and Pnth = 7200 K cm~^ 
(blue dashed curve) at the age of ~10^ yr (see the vertical dashed line in each panel). We also 
compare the SNR predictions with those for the LTC of our BPL (horizontal dotted lines in each 
panel with errors noted by the set of dotted lines at the top of right of the R12 and R45 panels; 
the LTC's O vi intensity and error bars are the same as those for the total BPL). The SNRs with 
no=0.02 atoms cm~^ (green curves) match those predicted by the LTC at age ~ 1.8 x 10^ yr (see 
the vertical dotted line), but the SNR models somewhat underestimate the LTC's R12 intensity 
(green solid and dotted lines) or overestimate its R45 intensity (green dashed line). 
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Table 1. IT and 2T Halo Model 





Local Bubble 


Halo (cool) 




Halo(hot) 






No. Model 


logT*'-'' EM''-'' 


logT''-'' EM''-'' 


rcf 


logT*-'' EM''-'' 




xVdof 



1 


IT CIE(S+R)'i 




8 9+°-'* 


6 •^5+0-01 












693.3/537 


2 


IT CIE(S)'' 


6.03 


8.9 




11.51°-, 










591.2/535 


3 


IT NEI(S) 


6.03 


8.3 


6 S7+0 01 


17 6+°-* 

^'•"-1.0 


35.O+3I 2 








573.8/534 


4 


2T CIE(S+R) 


5 92+0-04 


7 9+0.6 
'•^-0.4 


6 12+002 


24 l+''-'5 


6 t;n+o o2 


5 6+2-6 
^•"-0.7 




567.3/535 


5 


2T CIE(S) 


5.92 


7.2 


6 16+°°l 


I8.3I43 




6 ^1+0 03 
'J-'JJ^_0.02 


t; 1+0.8 




525.6/533 


6 


2T NEI+CIE(S) 


5.92 


7.2 




14.0l|3 


U.OUD_g 202 


6.51 


5.1 


1 •?6+*l*-® 
^-■^'^-0.76 


524.3/534 


7 


2T CIE+NEI(S) 


5.92 


7.2 


6.16 


18.3 




6 52+0-03 
"•^^-0.02 


8.01- 


519.4/534 


8 


2T NEI+NEI(S) 


5.92 


7.2 




"•''-2.6 


17+0-49 


fi '52+0 04 


7 9+2.0 
'•^-1.9 


01 7+18.3 
■^^■'-31.5 


516.7/531 



*In unit of K. 

''In unit of 10"^ cm"^ p^. 

'^In unit of 10^2 g cm~3, 

■^"(S+R)" means fit to Suzaku+Ft.OSAT data simultaneously. 

°"(S)" means fit to Suzaku speetra only. 

'The noted error bars reflect 90% confidence intervals. 



Table 2. Intrinsic Halo O Vll and O Vlll Intensities 



Model 




lo VIII 

(LU) 






9 c+O.S 



'^Obtained form the CIE 2T 
model fitting to the Suzaku spec- 
tra. 



'^Direct measurement, de- 
scribed in Section 5.1. 



- 33 - 



Table 3. Intrinsic Halo C iv, O vi, O vii, and O viii Intensities 



Ion 


Average Energy 


Intrinsic Intensity 




(eV) 


(LU) 



C IV 


-8.0 


7780 ± 2680 


VI 


-12 


7750t?g°o 


VII 


-570 


9.98^^^° 


VIII 


-650 


2.66t°j^ 



^The C IV and O vi error bars are la. The 
O VII and O viii error bars reflect the 90% confi- 
dence intervals. 



Table 4. Single Component P ower-law Halo DEM model Patterned on the Model in 



Shelton et al.1 (j2007|). 



logTz 


a 


VII 


VIII 


(K) 




(LU) 


(LU) 


6.06 


1.09 


2.28 


0.00 


6.24 


0.53 


10.4 


0.49 


6.36 


0.37 


18.4 


2.97 


6.54 


0.24 


22.9 


9.44 



Table 5. Fitting a Single Component Power-law Halo DEM Model to the Suzaku+ROSAT 

Spectra. 



Local Bubble 




Halo 










logT 


EM 


logTi 


logTa 


a 


VI 


-fciv 


xVdof 


(K) 


(10-3 cm-6 pc) 


(K) 


(K) 




(LU) 


(LU) 




o-yo_o.o6 

tr qr+0.05 

r q4+o.05 


6.4^^5 


5.06 
5.24 
5.57 
5.76 


6 fil+0 02 
6 fil+0 05 


-1 30+0-22 

-1 SO+0-21 

1 OO+0.27 
^■'J°-0.29 


46,800t^r 

46,500^ 
7850^1^^° 


83,000t^r 
18ltp 

38.0^23^2 


538.9/536 
538.9/536 
539.2/536 
541.9/536 



Note. — For each fit, logTi is fixed at the specified value, and the other halo and LB parameters are free to 
vary. The O vi and C IV intensities are then derived from the best-fit model parameters. The /dof pertains 
to the fit to the Suzaku+ROSAT data. 
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Table 6. Fitting Broken Power-law Halo DEM Models to the Suzaku+ROSAT Spectra. 





Local Bubble 


Halo{LTC) 


Halo(HTC) 






logT2 


logT 


EM 


logTi 




0-2 


logTg 




xVdof 


(K) 


(K) 


(10-3 cm-6 pc) 


(K) 




(K) 




(LU) 




5.97 
6.02 
6.06 
6.14 


5 03+0.04 

5 09+0 02 
■^•y^-0.05 
5 03+0.05 

c 00+0. 04 
O-»»_0.04 


fi s+0-5 
O-8-0.5 
fi s+0-5 

fi S+0-5 

fi 5+°-^ 


4.8 
4.8 
4.8 
4.8 


0.00 
0.30 
0.38 
0.54 


-2 10+°-" 

r, r,i +0.19 

^-^-^-0.12 

-2 06+°-^^ 

1 ^S+0-41 
--'--0»_o.l7 


6 

"■y^-0.02 
7 09+0 03 

6 96+° °'* 
O-»D_0.18 

"•y*-0.03 


9500 
7440 
6860 
5800 


544.2/536 
541.6/536 
541.8/536 
590.2/536 



Note. — For each fit, log Ti and log T2 are fixed at the specified value, and the other halo and LB parameters 
are free to vary. The C IV intensities are derived from the best-fit model parameters. The values of in 
final column are obtained by fitting the model to the Suzaku+ROSAT data. 



Table 7. Broken Power-law Halo Model: Testing Different Abundance Tables. 



Local Bubble Halo(LTC) Halo(HTC) 

Table logT EM logT\ logT^ ai logT^ /qvii J'oviii X^/dof 

(K) (10-3 cm-6 pc) (K) (K) (K) ( LU) ( LU) 



Grsa'' 


■s 05+0-05 


4 2+0-3 


4.8 


6.06 


0.11 


-2 16+°'i2 


fi Qfi+O OO 


11 1+^-'^ 


9 04+0. 63 


532.8/536 


Lodd'' 


K 05+0-05 


4-7ir3 


4.8 


6.06 


-0.25 


-2.40t»;ii 


fi qfi+0 00 
O^yD_0.02 


a fj7+2.60 
^•0'-1.63 


9 40+0-58 
^■^y-0.42 


541.1/536 


Wilm'^ 


5 QO+0.05 
''•='•-'-0.05 




4.8 


6.06 


0.38 


-2.06t»;?? 


6 qfi+0 04 


q fiQ+2-54 


2-47«;|I 


541.8/536 


Aneb<* 


c 02+0-04 
''•y^-0.05 


4 1+0-3 
^•^-0.3 


4.8 


6.06 


0.35 


-2 02+" " 


7 06+°'°° 


IO.8+1I 


2-86i°;Ii 


543.0/536 


Feld'= 


5 04+0 04 
"'•y^-0.05 


4 5+0-3 


4.8 


6.06 


0.72 




7 lfi+0-00 
'-^0-0.03 


10.6+?'« 


r,r,+0.63 

•5-'JU_o.44 


554.1/536 


Angr^ 


K 07+0-04 
•^•="-0.04 


4 0+0.3 
^••^-0.3 


4.8 


6.06 


0.60 


-2 05+'''l2 


7 Ofi+°'°° 


ll.lt?> 


9 01 +0.69 
^•°J--0.52 


563.3/536 



Note. — For each fit, logTi and logT2 are fixed at 4.8 and 6.06, respectively, and the other halo and LB parameters are free to 
vary. The O VII and O VIII intensities are then derived from the best-fit model parameters. The values of in the final column are 
obtained by fitting the model to the Suzaku+ROSAT data. 

HGrevesse k. Sauvall l ll998h 

iLodderd l l2003h 

' IWilms et al. I ll2000h . except that the XSPEC version sets several elemental abundance to 0. 

■ ^Anders fc Ebiharal l ll982h 

■ iFeldmanl lll992l') 

lAnders fc Grevessj lll989l ') 



Table 8. Broken Power-law Halo Model: Predicted Soft X-ray Count Rates 



Component R12(l/4 keV)'" R45(3/4 kcV)'' R67(1.5 keV) 





417+if 


-0.077 


U.UUOOO_Q ggQg4 


Low T Componenf^ 


339+1°' 


A 77+1-45 
' -0.86 


0.186lH^^ 


High T Componenf^ 


284l^f 


fjn f)+22.2 

' O.Z_;^3 2 


20.4l|;? 


EPL<= 


rr, 7+17.5 

^'■'-10.4 




104l?2 


Total 


10981333 


13314° 


124li 



''All values arc in unit of 10^^ ROSAT counts s^^ arcmin"^. 

'^Thc calculated SXR count rates for this component arc unabsorbed. 

'^The calculated SXR count rates for this component have been sub- 
jected to absorption due to TVh = 1-9 x 10^" cm~^. 



